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Abstract. 

We study chiral symmetry breaking in QCD, using as ingredients in the quark gap equation recent lattice results for the 
gluon and ghost propagators. The Ansatz employed for the quark- gluon vertex is purely non-Abelian, introducing a crucial 
dependence on the ghost dressing function and the quark-ghost scattering amplitude. The numerical impact of these quantities 
is considerable: the need to invoke confinement explicitly is avoided, and the dynamical quark masses generated are of the 
order of 300 MeV. In addition, the pion decay constant and the quark condensate are computed, and are found to be in good 
agreement with phenomenology. 
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One of the major challenges of the strong interactions 
is to understand the underlying mechanism that generates 
masses for the quarks and breaks the chiral symmetry 
(CS). The CS breaking is an inherently nonperturbative 
phenomena, whose study in the continuum leads almost 
invariably to a treatment based on the Schwinger-Dyson 
(SD) equation for the quark propagator (gap equation). 
As is well known to the SD experts, the existence or not 
of solutions for this equation depends crucially on the 
details of its kernel, which is essentially composed by 
the fully dressed gluon propagator and the quark gluon 
vertex [1]. The latter quantity controls the way that ghost 
sector enters into the gap equation, and introduces a 
numerically crucial dependence on the ghost dressing 
function [2] and the quark-ghost scattering amplitude. 

In the present talk we report on a recent study of CS 
breaking [1] using the SD equation for the quark prop- 
agator, supplemented with three nonpertubative ingredi- 
ents: (i) gluon propagator and (ii) ghost dressing function 
obtained from large-volume lattice simulations, and (iii) 
the "one-loop dressed" approximate version of the scalar 
form factor of the quark-ghost scattering kernel. 

The starting point is to express the fully dressed quark 
propagator in the following general form [3] 



S- 1 (p) => -m - Z(p) =A(p 2 )rf -B(p 2 
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FIGURE 1. The quark SD equation (gap equation). 



The SDE for the quark propagator, which is repre- 
sented schematically in Fig. 1, can be written as 



S- l (p) =j/-C F g z Jj^S(k)T v (-p,k,q)^ v (q), (2) 

where q = p — k, J k = pL 2£ (2%)- d Jd d k, with d = 4 - e 
the dimension of space-time. Cp is the Casimir eigen- 
value in the fundamental representation (for ST/ (3) 
Cp = 4/3). The full gluon propagator A liv (q), in the Lan- 
dau gauge, has the form 



A" v ( 9 ) = -i 
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A(<7 2 ), 
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where the non-perturbative behavior of the scalar fac- 
tor A(q 2 ) has been studied in great detail in the con- 
tinuum [4, 5] and lattice simulations [6, 7]. The fully- 
dressed quark-gluon vertex T v (—p,k,q) also obeys its 



where m is the bare current quark mass, and T.(p) the 
quark self-energy. We consider the case without explicit 
CS breaking, i.e., bare mass m = 0. The dynamical 
quark mass function can be defined as being the ratio 
^£(p 2 ) =B(p 2 )/A(p 2 ). Then CS breaking takes place 
when the scalar component, B(p 2 ) develops a non-zero 
value. 




H(puP2,P3) = P 1 
FIGURE 2. The quark-ghost scattering kernel. 




FIGURE 3. Lattice results for the gluon propagator, A(q), and ghost dressing, F(q) renormalized at p. = 4.3 GeV. 



own SD equation which unfortunately, it is too compli- 
cated. Then, we have to resort to the so called gauge- 
technique, where a nonperturbative Ansatz for the vertex 
F ll (pi,p2,p3) written in terms of S(p) is constructed, 
based on the requirement that it should satisfy the fun- 
damental Slavnov-Taylor identity (STI) 



pfrp = F( P3 ) [S- 1 ( -p , )H - HS 1 fa)] 



(4) 



where the ghost dressing function F(p 3 ) is related to the 
full ghost propagator D(p 3 ) by D(p 3 ) = iF {]}■$)/ p^. The 
quark-ghost scattering kernel H, represented in Fig. 2, 
and its "conjugated" H are functions of the momenta 
H = H(pi , p 2 , P3), H = H_(p 2 ,pi,pi) respectively. 

Both kernels H and H have the following Lorentz 
decomposition [8] 

H{pi,p 2 ,p3 hX I + Xrf l+X2jfc+X 3 OnvrfPl , 
H(p 2 ,Pi,P3)=X I-X^ -X^ 2 +X 3 a^pl (5) 

where the form factors X, are functions of the mo- 
menta, X, ■ = Xi(p\,p2,pj), and we use the notation 
Xi{p2,P\,P3) =Xi(p\,p2,P3) and a^ v = jh^/v]. 

On the other hand, the most general Lorentz structure 
the longitudinal part of the vertex F ll (pi,p2,pi) appear- 
ing in the lhs of Eq. (4) is given by [8] 

^fi{Pl,P2,P3) =L 1 Yl 1 +L 2 (p' 1 -fa)(pi-pz)p 

+L i (pi-p 2 )^+L4d llv (pi-p2) v , (6) 

where once again we have suppressed the depen- 
dence on the momenta in the form factor L, [i.e. 
L{ = Li(p\ ,p2,P3)] -Notice that, the tree level expression 
is recovered setting L\ = 1 and L,2=L 3 =L\ = 0. 

Due to the fact that the behavior of the vertex is 
constrained by the STI of Eq. (4), the form factors L;'s 
appearing into the Eq. (6) will be given in terms of the 
form factors X, 's of Eq. (5). 

The full expressions for L,'s in terms of the form 
factors Xj's is given in [1]. For the sake of simplicity, we 



will show here the case where only the scalar component 
of the quark-ghost scattering kernel is non-vanishing i.e. 
X ^ while Xi =X; = 0, for i < 1. In this limit, we 
obtain the following expressions 
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According the above expression the form factor 
Li's displays an explicit dependence on the product 
F (p 3 )Xo(p 3 ) which contains information about the IR 
behavior of the ghost propagator. Therefore, the ghost 
sector couples to the gap equation Eq. (2) through the 
quark-gluon vertex of Eq. (6). It is interesting to notice 
that in the limit of F(p 3 ) = Xo(p 3 ) = 1 the form factors 
of Eq. (7) reduces to the ones used so-called Ball-Chiu 
(BC) vertex [9]. 

We will next insert into Eq. (2) the general quark- 
gluon vertex of Eq. (6) with the expressions for the form 
factors Lj given in Eq. (7). Defining p\ — —p, p2 = k, 
and p 3 = q and taking appropriate traces, we derive the 
expressions for the integral equations satisfied by A{p 2 ) 
and B{p 2 ) that schematically can be written as 



A(p') = l+C F g%. 
B( P 2 ) = C F g 2 Z- 1 



■Mp-k) 



k A 2 (k 2 )k 2 +B 2 {k 2 ) 
Jto(p-k) 



Jt A (k,p), 



k A 2 (k 2 )k 2 +B 2 (k 2 ) 



^b{Kp), (8) 



where the kernel j£o(q) corresponds to the part that is 
not altered by the tensorial structure of the quark-gluon 
vertex, namely J^o(q) — A(q)F(q)Xo(q), while the parts 
affected are J^(fc,/?) and Jtg(k,p). 

The gap equation depends on the nonperturbative form 
of the three basic Green's functions, namely A(q), F(q), 
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FIGURE 4. The form factor X^ l] (q) given by Eq. (9 



and Xo(q). For A(q) and F(q) we use the recent lattice 
data obtained by [6], and shown in Fig. 3. 

We clearly see that both lattice results for A{q) and 
F(q) are infrared finite. Such a feature can be associated 
to a purely non-perturbative effect that gives rise to a 
dynamical gluon mass [10], which saturates the gluon 
propagator in the IR. The appearance of the gluon mass 
is also responsible for the infrared finiteness of the ghost 
dressing function, F(q 2 ) [4, 11], which is shown on the 
right panel of Fig. 3, 

Unfortunately forXo(g) there is no lattice data avail- 
able, and in order to obtain a non-perturbative estimate 
for Xq, we will study "one-loop dressed" scalar con- 
tribution of the diagram of Fig. 2 in an approximate 
kinematic configuration, which simplifies the resulting 
structures considerably. Specifically, we will assume that 
pi = p2 = p, and p — —q/2. In doing so, we arrive at 
(see details in [1]) 



X, 



1 



(q) = l + -C A g 2 q 2 



(k-qf 



A{k)F{k) 



F[k + q) 

(k 



f<7) 4 
(9) 

We proceed substituting the fit for the lattice data for 
A(q) and F(q) presented in Fig. 3 into Eq. (9). The 

numerical result for xjp (q) is shown in the Fig. 4. 



Dynamical quark r 

BC vertex 

Asymptotic behavior 




X^\q) shows a maximum located in the intermedi- 
ate momentum region (around 450 MeV), while in the 
UV and IR regions X^\q) — > 1. Although this peak is 
not very pronounced, it is essential for providing to the 
kernel of the gap equation the enhancement required for 
the generation of phenomenologically compatible con- 
stituent quark masses. 

Now we are in position to solve the sys- 
tem formed by Eq.(8) Substituting A(q 2 ), F(q 2 ), 



and 
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P[GeV] 



to Eq.(8), with the modification 
Z~ L j£ AB (k,p)^J(f A . B (k,p)F(p 2 ), to enforce the 
correct renormalization group behavior of the dynamical 
mass (see discussion in [1]), we determine numerically 
the unknown functions A(p 2 ) and B{p 2 ). The result for 
the dynamical quark mass ^{p 2 ) is shown in Fig. 5. 

One clearly sees that ^{p 2 ) freezes out and acquires 
a finite value in the IR, J?(0) = 294 MeV. In the UV it 
shows the expected perturbative behavior represented by 
the blue dashed curve. 

With the behavior of the dynamical quark mass at 
hand, we have computed pion decay constant and the 
quark condensate and we obtained f K = 80.6 MeV and 
(qq) (IGeV 2 ) = (217MeV) 3 respectively, which are in 
good agreement with phenomenological results. 
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FIGURE 5. The quark mass . 



i for the BC vertex. 



